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Recently discovered models of ghost-free massive gravity and bigravity are characterized by a 
non-trivial potential that gives rise to a rich vacuum structure. We review maximally symmetric 
vacua of the de Rham-Gabadadze-Tolley (dRGT) massive gravity and of the Hassan-Rosen (HR) 
bigravity, and discuss their perturbative stability. In particular, we discuss perturbations about 
self-accelerating vacua in HR bigravity, and argue that, analogously to what was found in the case 
of dRGT gravity, some of them contain strongly coupled modes. We then show that it is impossible 
to construct regular instantons connecting different classically stable vacua of dRGT gravity without 
violating energy conservation or the null energy condition. 
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I. INTRODUCTION 

A scries of papers [1-7] have recently led to the construction of models of massive gravity (denoted as "dRGT" 
from the authors of [8]) and bigravity ("HR" [9]) that propagate only five (or seven, in the case of bigravity) degrees 
of freedom and are therefore free of the Boulware-Deser (BD) ghost [10]. Remarkably, these models are described 
by a finite (ranging from 3 to 6, depending on the assumptions) number of parameters. Such a parameter space is 
large enough to yield a rich vacuum structure, while being small enough to allow for a complete characterization of 
the vacua of the theory. Besides the obvious Minkowski vacua, various cosmological solutions [11-19] and spherically 
symmetric backgrounds [20-23] have been discovered. 

In the present paper, complementing previous works, we characterize the landscape of the SO(4)-invariant vacua 
of Euclidean dRGT massive gravity and investigate their perturbative and - to our knowledge, for the first time 
- non-perturbative stability. As we will see, we will recover two disjoint branches of solutions: vacua in the form 
of constant factor x Minkowski and (Anti-)de Sitter-like vacua. Perturbative analysis shows that depending on the 
choice of parameters some of these vacua will be stable, while others will be plagued by tachyonic and/or ghost-like 
pathologies. 

Given the presence of multiple vacua, it is natural to ask whether there can be non-perturbative instanton transitions 
bringing one vacuum to a different one. We will show that it is impossible to construct solutions describing such 
a transition unless we allow for a non-conserved source or for the violation of the null energy condition in the 
matter sector. These results are similar (but, as we will see, the details are quite different) to the findings of [24], 
who studied the possibility of non-perturbative transitions in the context of the Dvali-Gabadadze-Porrati (DGP) 
model [25]. Zhang et al. [26] also looked into instantons in the context of dRGT gravity. However, these authors 
studied instantons between the different vacua of a given scalar field, whereas in the present work we are concerned 
with the vacuum structure of dRGT gravity itself. It is worth pointing out that the quantum consistency of the 
models under consideration has not been fully explored yet, and there might be inconsistencies. For example, the 
structure found in the papers [8, 9] that guarantees the absence of the BD ghost may not be stable under radiative 
corrections [27, 28]. Moreover, like any theory of massive gravity, these models suffer from strong coupling at low 
energies that would cause a loss of predictability at scales as large as a kilometer [29] . Since instanton transitions are 
a purely quantum effect, our effort should be understood as another probe to the quantum aspects of massive gravity. 

For completeness, we will also extend our program of vacuum search to HR bigravity, finding again two classes of 
maximally symmetric vacua, which in the limit where the second graviton decouples converge to the corresponding 
vacua of dRGT. Then at perturbative level, we will study the dynamics of vector modes on one such class of vacua: 
self-accelerating vacua in dRGT are known to contain strongly coupled modes [30-33]. Since dRGT corresponds to 
HR in the limit where the second metric becomes non-dynamical, one might ask whether strong coupling is an artifact 
of the decoupling of the second metric. It turns out not to be the case: even when both metrics are dynamical, the 
self-accelerating branch contains (infinitely) strongly coupled vector modes, consistently with the findings of [34]. We 



*minjoonp@physics .umass . edu 
tsorboOphysics.umass.edu 



2 



also study the stability of perturbations about vacua of the form of constant factor x de Sitter in HR gravity. We 
leave the investigation of instantons in HR to future work. 

The plan of the paper is the following. After introducing a general formulation of massive (bi)gravity theories in 
the next section, in §111 we will find the S'0(4)-symmetric Euclidean vacua of dRGT, and calculate the linearized 
action about them to obtain criteria for perturbative stability. Then in Section IV, we check the non-perturbative 
stability of dRGT by trying to construct instantons for all possible configurations of vacuum transitions. Section V 
is devoted to the background and the perturbative analyses of HR theory, and we conclude in §VI. Heavy algebraic 
details can be found in the Appendix. 



II. FORMULATION OF GHOST-FREE MASSIVE (BI) GRAVITY 

The most general form of BD ghost-free massive gravity is the Hassan- Rosen (HR) bigravity [9] , where two metrics 
g^v and are both dynamical and whose action reads 



(1) 



S=^J d 4 xy^g- W R +^Lj d 4 x^f if) R + m 2 M 2 s J d A x^~g L + S, 



where M ofi 2 = M~ 2 + MJ 2 , 7 £ 7 £ = g^f pv and 

4 



n=0 



+ _i [ 7 4 _ 6 7 2 7 . 7 _|_ 877 • 7 • 7 + 3(7 • i) 2 - 67 • 7 • 7 • 7] 



(2) 



In order not to reintroduce the BD ghost, matter should couple to either one of the metrics but not to both, and we 
assume only g couples to matter. Using 



we get the full EOMs: 
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with T£ the stress-energy tensor for matter and 



T* = ]T Pne [n f v {l) = M + fcillZ ~ iW) + f frV 277p^7^ + 



n=l 



+ § [7(7 2 ~ 3 7 • 7)7^ + 3(7 • 7 - 7 2 )7^7" + 2 7 • 7 • 77^ + 677^7^ - ^Ipl^. 



f v = ^Pi-ne^il 1 ), Lp = ^/3 4 -„e ( „)(7 x ) . 



(6) 



71=1 



n=0 



To obtain the dRGT model, we first decouple the reference metric by taking Mf — > 00, which also results in 
M c ff — > M g . In this regime the term in /?4 does not contribute to EOMs and (5) is solved by any Ricci-flat f^. For 
generic values of (3 , • • ■ , P3, however, the equations of motion do not allow for a Minkowski solution g^ u = r]^ unless 
Po + 3 (3i + 3 (3 2 + /?3 = 0. In this case, it is customary to redefine the parameters as 



6-4a 3 +a 4 , ft = -3 + 3a 3 -a 4 , (3 2 = 1 - 2 a 3 + a A , /3 3 = a 3 - a 4 . 



(7) 



Thus the dRGT model is defined with three parameters, a 3 , 0:4 and to, where a's give various interactions between 
/ and g, and to is a parameter with the dimension of mass that will turn out to give the mass of the graviton about 
the vacuum with g^ v — -q^. 



III. DRGT GRAVITY 
A. Vacua 

We study the vacua of the dRGT theory in Euclidean signature and we confine ourselves to SO(4)-symmetric 
metrics. Our Ansatz is then 

g^dx^dx" = a{rf dr 2 + b{r f dfl 2 u , (8) 
f^dx»dx v = dr 2 +r 2 dn 2 m , (9) 

where we assume a, b > 0, and dil 2 u = d\ 2 + sin 2 \ df2 2 j = d% 2 + sin 2 x(d^ 2 + sin 2 9 dip 2 ) is the metric on S 3 . Using the 
Ansatz (8-9), the equation of motion (4) has only two independent components. But one of them being redundant, 
we have only 

3bb' 2 -3ba 2 - m 2 [(6 -4a 3 + a 4 )b 3 - 3 (3 - 3 a 3 + a 4 ) r b 2 + 3 (1 - 2 a 3 + a 4 )r 2 b+(a 3 - a 4 )r 3 ] a 2 = . (10) 

The other equation comes from the Bianchi identity, i.e., the covariant divergence of (4) should vanish: 

(a -b') [(3- 3a 3 + a 4 )b 2 - 2 (1 - 2 a 3 + a 4 ) r b - (a 3 - a 4 ) r 2 ] = . (11) 

Eq. (11) can be solved by either a = b' or b = x±r, where 

_ 1 - 2a 3 + a 4 ± y/1 - a 3 + a\ - q 4 
6 — 6a 3 + a 4 

and we have the following two branches of solutions: 

1. Cosmological Solution (CS): For b = \± r i solving (10) for a gives 

d 

1 - fi{x±)m 2 r 2 ' 

with 



1 — x 

H{x) = — — [(6 - Aa 3 + a 4 ) x 2 - (3 - 5 a 3 + 2 a 4 ) x - a 3 + a 4 ] , (14) 



and the physical metric reads 



g^dx^dx" = - 2 dr 2 + X \ r 2 dfijf n . (15) 



1 - Mx±) ™ 2 r 

With a double Wick rotation we write the metric in Lorentzian signature as 

Y 2 

g^dxW = -- X± „ ? dr 2 + X 2 ±T 2 (dR 2 + sinh 2 i?dtt 2 T ) , (16) 

1 + fi±m z t 

with /i± = /Li(x±)- Since /i + > 0, we can bring the metric with \+ to a cosmological form through a redefinition 
r 



"»V7*+ x+ 

^dxMx" = -dt 2 + ( X+ sinh TO V^+ A 2 ( dfl 2 + sinll 2 fld Q2 J _ 

This is the de Sitter metric in open chart found in [11]. As for the solution with x_, since ^_ < 0, it turns out 
to describe an Anti-dc Sitter metric* 



'Arbitrary cosmological solutions were found in [15, 16], where took complicated forms. The de Sitter solutions found in those papers 
would agree with the open solution of [11] up to coordinate transformations. For example, after a standard coordinate transformation 
from open dS to flat dS, the vacuum solution of [11] can be matched to the de Sitter solution of [15] with /(t, r) = ^-Hr 2 e Ht + ^ sinh Ht, 

9 2 

where H z = Pq(xq), and Pq> / and xq are defined in §111 of [15]. 
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2. Conformal-to-Minkowski (CM): If a = b' , g^ v is fiat, and eq. (10) has three solutions: b = r and b 

_ 3 - 5 a 3 + 2a 4 ± ^9 - 6a 3 + 9a§ - 12a 4 
C± ~ 2(6 - 4a 3 + a* 4 ) ' 

Therefore we have a Minkowski and two conformal-to-Minkowski vacua: 

g^dx»dx v = c 2 (dr 2 + r 2 d0 2 n ) , (18) 

where i = +, 0, — and cq = 1. For all three solutions, the physical metric g^j, can be transformed into Minkowski 
with a simple coordinate rescaling. But of course such a rescaling causes the fiducial metric to move away 
from Minkowski, so that each solution corresponds to a different physical situation; while ordinary matter will 
behave exactly the same on the three backgrounds with b = c+ r, b = c_ r, and b = r, graviton fluctuations will 
feel a different mass. For example, as we will see below, when 6 — 4 a 3 + a 4 > 0, perturbations about the metrics 
with Co or c_ describe a graviton with positive mass, whereas perturbations about c + -vacua are tachyonic. 



B. Perturbations of the CM backgrounds 

Let us now perform a perturbative analysis on top of CM backgrounds described above and investigate their 
stability.''' Since only g is dynamical, we expand (1) as g^ v — cf rj^ + h^, where x^ = (t, x) and 

h i = hl + diB, htj = hJ? + d i Sj + d j &-26 ij il; + 28 i d j E, (19) 

with difiQi = 0, SijhJ^ = 0, dihj? — and d&J — 0. To handle the expansion of 7 = g~ x f, we use the trick of [31]: 
For N x N matrix A, with 

A = diag(oi,--- ,a N ) + eA^ + e 2 A™ +0(e 3 ), a t > , (20) 

VA = diag(^I, • • • , ^/on) + eB« + e 2 B< 2 ) + 0(e 3 ) , (21) 

B (i). = J_ j B (2), = V ^ _ (22) 



its square-root matrix 



is given by 



Since different helicities do not mix at the quadratic level of the expansion of (1), we consider each helicity mode 
separately. The quadratic action for helicity-2 modes is easy to obtain, whereas for helicity-1 and -0 modes, we have 
to integrate out/solve for non-dynamical modes (equivalently, one can take Mf — > 00 limit of the derivation of the 
quadratic action of HR bigravity presented in the Appendix): 

Shd- 2 = ^ / d* 4 {\(W) 2 + 4^ T ( A - -cW} > ( 23 ) 

= f / ^{\^^Jf - \^i^Jf} , (24) 

Shd-o = ^-J dx 4 {6^ 2 + 6^(A-m 2 ff )^}, (25) 
where A = &f and where 

2 / n 2 /l-2a 3 + a 4 2(3-3a 3 + a 4 ) \ . . 

mi s {ci) = m A y -2 - h 6 - 4a 3 + a 4 J . (26) 



It is now obvious that the sign of m 2 ff determines the stability of the model. By plotting m 2 ff for a 3 and a 4 , one can 
see that 



t Perturbations in the CS branch were investigated in [31, 32, 35, 36]. 
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FIG. 1: Marked by blue shade is the region in the a-plane where two stable CM vacua are allowed. The two straight lines are 
tangent to the parabola. 



• TOg ff (l) = 1, so that a theory with a Minkowski vacuum is well behaved for any value of 0:3 and a 4 , at least at 
the linear level. 

• m^ ff (c+) is negative, and therefore in a model built on the c+-vacuum, all the modes, (23-25), are tachyonic. 
Furthermore, the helicity-1 mode becomes a ghost at short enough scales. 

• When 6-4«3 +«4 < 0, even m^g(c_) becomes negative, so that c_-vacuum as well as c+-one has tachyonic 
instabilities. 

These results can also be explained by the following argument: If we look for metrics conformal to Minkowski, 
9fiv = < f 2 Vfiui we see that <E> has a quartic potential 

V($) oc $ [(6 - 4 a 3 + a 4 ) $ 3 - 4 (3 - 3 a 3 + a 4 ) $ 2 + 6 (1 - 2 a 3 + or 4 ) $ + 4 (a 3 - a 4 )] , (27) 

whose extrema are located at $ = c + , I and c_. For 6 — Aa 3 +a 4 > 0, it can easily be seen that e_ < c+ < 1, so that 
$ = I and c_ are local minima, while $ = c + is a local maximum and therefore unstable. On the other hand, when 
6 — 4 a 3 + a 4 < 0, the ordering among CM vacuum solutions changes into c + < I < c_ . Combining this with the fact 
that the potential (27) is now unbounded from below, we see that the Minkowski vacuum is the local minimum and 
c± -vacua are local maxima, hence the appearance of instabilities for both c+ and c_ backgrounds. 

From the perturbative analysis, we may conclude that if we want dRGT to allow for two stable CM backgrounds 
(c_ as well as Co), we should constrain the parameter space of (03, a 4 ), such that 

6 - 4a 3 + a 4 > 0. (28) 

We also assumed b > at the beginning, which requires c± be real and positive. Under the restriction (28), this 
implies 

3-2a 3 +3a3-4o!4>0 and a 4 - a 3 > . (29) 

Therefore, while Minkowski is always a stable solution, the second, stable, conformal-to-Minkowski vacuum is allowed 
only on the intersection of (28) and (29), which is depicted in Fig. 1. 
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IV. INSTANTONS IN DRGT 



Since dRGT gravity can have multiple classically stable vacua, it is natural to ask whether non-perturbative 
solutions of the equations of motion can allow for transitions between them. As long as we limit ourselves to the 
SO(4) symmetric solutions of pure massive gravity without matter source, it is clear that no such transitions can 
occur - the metrics found in §111 A include all possible solutions for such a system. However, matter might be able to 
act as a catalyst for vacuum transitions. In this section we will first discuss the possibility that some singular form of 
matter supports the instanton. We will also show that, as long as matter satisfies the null energy condition and the 
energy momentum tensor is conserved, there are no smooth solutions describing transitions between the various vacua 
described above. Given that we are considering vacua belonging to separate branches, and that we will be discussing 
instantons supported either by singular or by regular matter source, this section is divided into several subsections. 



A. (No) Singular instantons within the CM branch 



A singular metric interpolating between two vacua in the CM branch should be written in a form that allows 
the equation of motion (4) to make sense at least in a distributional form. As a consequence, the presence of the 
Einstein tensor built out of g^ v requires this metric to be continuous. Since we are considering transitions within 
the CM branch, we must therefore choose the gauge where g^ v is exactly Minkowski both "inside" and "outside" the 
instanton, whereas we allow the auxiliary metric f^ v to have some form of discontinuity across the instanton. As a 
consequence, we might look for solutions of the form 

6V = dr2 + r2 d0 ni i 

_ f cZ 2 (dr 2 + r 2 dn 2 m ) , r<r Q , . 

■>^~\dr 2 + r 2 dn 2 m , r>r ' [M) 

(We are not interested in any transition involving the c + -vacuum, because it is classically unstable.) However, 
the metric / M „ in eq. (30) carries curvature on the singular surface r — r and does not correspond to dRGT 
gravity that requires f^ v to be identically Ricci-flat. Then, one can try f^dx^dx 1 ' = y'(r) 2 dr 2 + y(r) 2 dil 2 u with 
y(r) = cZ 1 r + (1 — cZ 1 ) r 6(r — r ). Such f^ v is flat everywhere by construction, but it contains a singularity 
cx 5(r — r ) 2 in f rr , inducing a term proportional to S(r — r ) in 7^. While this by itself is not a problem, since 
the theory contains powers of 7^, there will be terms with powers of the Dirac delta, which are not meaningful in a 
distributional sense. Equivalently, one might go back to the Stiickclbcrg formulation of theory by introducing a set 
of four scalar fields <f) a , such that f^ v — <5 ab <9 M </>° d v <j> h , and 





= K 




cZ^Oir- 


ro)} 


r cos x , 


? 


= K 




cZ^Oir- 


ro)} 


r cos 9 sin x , 


? 


= K 


'+(1- 


cZ^Sir- 


ro)} 


r sin 9 sin ip sin \ 


, 4 


= K 




cZ^Qir- 


ro)} 


r sin 9 cos f sin x 



where x, 9 and if are the angles on S 3 . Again, since functions <j) a are discontinuous, powers of the Dirac delta function 
will appear in the equations of motion (4) , making them ill-defined from a distributional point of view. 

Therefore we conclude that the theory does not allow for singular and yet meaningful in a distributional sense, 
non-perturbative solutions describing transitions between different CM vacua. 



B. (No) Singular or regular instantons within the CS branch 



In the CS branch a given choice of a 3 and 0:4 may allow two different vacua; a de Sitter vacuum with b = x+ r and 
an Anti-de Sitter one with b — X- r - An instanton connecting x + -vacuum to x--one looks like 

g^dx^dx" = ai„(r) 2 dr 2 + X in(r) 2 r 2 df> 2 „ , (31) 

where a In jumps from , x+ = to , x ~ ==, while Xin goes from x+ to X-- Now, as in the previous case, we 

have to choose a coordinate system where the induced metric at the junction of the two geometries is continuous, 
but in doing so we generate a singularity in the / M „ metric that makes the equations of motion ill-defined from a 
distributional point of view. 
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We can go further, and prove that no regular transition, supported by matter that obeys energy conservation, can 
happen within the CS branch: the Bianchi identity (11) requires b should always be either \+ r or X- r , *- e -> Xin be 
cither \+ or X-- If wc require a smooth metric, then \in cannot change at all. Note that the form of the matter does 
not play any role (as long as the stress-energy tensor is conserved) in this argument, and there can be no singular or 
regular instanton connecting different vacua within CS branch. 



C. (No) Singular transitions between vacua of different branches 

For the case of a transition between a vacuum in the CM and one in the CS branch, we can use an argument 
identical to that of section IV B above. Let us consider for instance the case where we are in the CM branch, 

9flv dx»dx» = c 2 (dr 2 + r 2 dSl 2 u ), f^dx^dx" = drf + r 2 dft ni , n < f< , (32) 

at small radii (r^ being a radial coordinate covering the "interior" of the instanton, and Cj denoting one between 
c (= 1) and c_), whereas at large r we are in the CS branch, 

g M „ d/d/ = _ X " 2 2 dr 2 + xl r 2 a dft nl , /„„ dx»dx v = dr 2 + r 2 dft nl , r > f , (33) 

where r covers the "exterior" of the instanton geometry, and \n is either \+ or X-- 

Then by requiring the metric g^ to be continuous at the junction of the two geometries we must impose fj = Xn^o, 
which, however, implies a discontinuity f 2 (l — cf/x 2 ) in the coefficient of the df2 2 n term in the non-dynamical metric 
f^. It is easy to see (see also the discussion in the next subsection) that in general 7^ Xn- Therefore we are forced 
to introduce a discontinuity in / M „ that makes the equations of motion ill-defined in a distributional sense, and hence 
singular transitions between vacua of different branches are not allowed. 



D. (No) Regular vacuum decay within CM branch 

Since the various vacua of the CM branch all belong to the same solution, a = b', of the constraint (11), it is 
possible to transition among them without violating energy conservation. Therefore we have to analyze in detail the 
dynamical equations to see what conditions must be satisfied by the matter supporting the instanton. 

An Ansatz for an instanton within the CM branch can be written as 



g^dx^dx" = 6j n (r) 2 dr 2 + 6 In (r) 2 d^ n; (34) 
f^dx^dx" = dr 2 + r 2 dn 2 u . (35) 

Then (4) gives 

T r p m 2 

-TJ2 = ~T/2 = I5~ {r-hn) [(6 -4a 3 + a 4 )b 2 n - (3 - 5a 3 + 2a 4 )r6 In - (a 3 - a A )r 2 } , (36) 

9 9 In 



Ml = Ml = blb> ln 



— < [3 - 3 a 3 + a 4 - (6 - 4 a 3 + a 4 ) b' In ] bj n 

'in L 

-2[l-2a 3 + a 4 - (3 - 3a 3 + a 4 ) b' In ] rb ln - [a 3 - a 4 + (1 - 2 a 3 + a 4 ) b' ln ] r 2 } , (37) 



where p and p are the energy density and the pressure of the matter supporting the instanton. Defining c(r) = b\ n (r) /r, 
we get 

2 f 

-^ff = "^sTT- [(3-3a 3 + a4)c 2 -2(l-2a3 + a4)c-a3+a 4 ] , (38) 

M g c "in 

where the left hand side has to be positive in order to satisfy the null energy condition. * Note that the zeros of 
s(c) = (3 — 3 a 3 + a 4 ) c 2 — 2 (1 — 2 a 3 + a 4 ) c — a 3 + a 4 are at c = x± j an d one can easily show that within the allowed 



^For a canonically normalized scalar field <j> with potential V(<j>), wc would have p = — <ji' 2 /2 6j n 2 + V(4>) and p = — tf>' 2 / 2 b' ln 2 — V(<f>), so 
that the LHS of (38) equals (j>' 2 / M 2 b[ n 2 > 0. 
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region, (28-29), of a's, 

c_ < x- < c+ < x+ < 1 , (39) 

i.e., 

s(c) { > ' for c < X- or c> X + , m 

S[C) \ <0, for X - <c< X+ . [W) 

We now have all the necessary elements to show that it is impossible to find a healthy solution where c(r) interpolates 
between c_ and 1 as r ranges from to oo. 

Let us first assume that c(r — > 0) — > 1. Then since s(l) = 1, it follows that c'(r — > 0) > 0, i.e., if c starts at 1, it 
can only increase. Therefore c will never be able to decrease to reach c_ < 1. This excludes the possibility that the 
interior of the instanton is the vacuum with c = 1. 

When c(r — > 0) — > c_, in order to reach c = 1, c(r) needs to cross first X- an d then x+- If we do not want to 
violate the null energy condition, the quantity c' (r)s(c)/b' ln (r) should stay positive as c crosses, say, X-- Now, c(r) 
might have wiggles, and can cross \- many times, but let us define r_ by the smallest among r's solving c(r) = x_. 
Since, at r = r_, c(r) crosses X- from below, c'(r_) > and therefore &j n (r_) = c'(r_)r_ + c(r_) > 0. Then c'/b' ln is 
positive at least in a small neighborhood of r = r_. But in the same neighborhood, s(c) changes signs from positive 
to negative as c crosses X- ; and so does c' (r)s(c) /b' ln (r), violating the null energy condition. 

Therefore we see that there is no regular instanton that can connect the two stable CM vacua. 

E. (No) Regular transitions between vacua of different branches 

We would like to check whether it is possible to smoothly connect a CS background, 

g^Wdx" = ^^dr 2 +xlr 2 <m 2 m , n = ±, (41) 

at, say, small r, to a CM one, 

g^dx^dx" = c\ (dr 2 + r 2 d0 2 „) , i = 0, - , (42) 
at large r. Since generic CS and CM solutions in the presence of a source can be written as 

g^&x^&x" = a(r) 2 dr 2 + xl r 2 d^ 2 „ and g^dx^dx" = b'(r) 2 dr 2 + b{rf dft 2 „ , (43) 
respectively, the solutions we are looking for should behave as follows: 

2 

• As r increases from to a certain r*, a{r) 2 smoothly deforms from 1 _^ x ^ 1 2 r 2 to x 2 > so that at we reach 
g^dx^dx" = Xndi" 2 + Xn 7,2 which is both CS and CM. That is, at r — r*, we can also say we have a CM 
solution with b — Xn T - 

• Then, as r grows from r» to 00, b(r) smoothly transforms from Xn r to Cir. 

With four vacua (x+, X-; c o an( l c_) to decay to and from, we need to consider 8 different cases: 



~— — - — r > 
r < r* 


co 


c_ 




X- 


co 




I 


II 


c_ 


III 


IV 


X+ 


V 


VI 




X- 


VII 


VIII 



Before we investigate each of them, let us list some features of CM solutions, which will be useful for our investigation. 

1. From the analysis of §IVD, we saw that if c(r — > 0) — > c (= 1), the null energy condition requires c'(r — > 0) > 
and c monotonically increased. 



9 



2. Similarly, one can show that if c(r — > oo) — > c_, c'(r — > oo) should be positive in order to satisfy the null energy 



3. In §IVD, we also observed that when trying to connect c_ and 1 (or even x+) with such c(r) that c(r — > 0) — > c_, 
the null energy condition was violated as c crossed X- ■ 

4. Similar violation of the null energy condition will occur when c crosses x+> if we try to connect \- an d 1 with 
c(r) satisfying c(r — > oo) — > 1. 

Now it is easy to see that there is no instanton solution for I and II: For r < r t , we have a CM solution with 
b(r) = c(r)r and c(r — > 0) — > 1. Then by item 1 of the list above, c keeps increasing as r increases, being unable 
to reach \± (< 1) at r = r*. Likewise, item 2 forbids any transition of type VI and VIII, because c(r) would only 
decease from c_ (< x±) as r deceases from oo. It is also obvious that no vacuum decay of type III and VII can occur, 
because of item 3 and 4, respectively. 

To examine IV, we first obtain that for a CS solution, — (p + p) oc Xn ra ' + Xn a — a3 > so that NEC demands 
X^ra' > a(a 2 — x«)- That is, if a(r*) = Xn, a(r) will keep increasing (deceasing) as r gets larger and larger (smaller 

2 

and smaller) than r*. For the case of IV to be realized, a(r) 2 should interpolate between x~ at r = and 1 _J C m 2 r 2 

2 

at larger r. But since /u_ < 0, 1 _J C m 2 T -i < X- f° r an Y r ; whereas a(r) gets bigger than X- as r grows, so that a(r) 
cannot be deformed into m 2 r 2 , making it impossible to have an instanton of type IV. Similarly V cannot be 

2 

realized either: 1 _^ + + t2 r 2 > x\ because /x + > 0. But as r deceases from r*, a(r) gets smaller than x+- 

Therefore, we can conclude that there is no regular vacuum transition between vacua from different branches. 



The discussion in §IVD concerns transitions between two different stable CM vacua. For certain values of the 
parameters «3 and a.±, however, there is only one possible vacuum. This may occur either 

• when 6-4a 3 + a 4 > and 3 — 2 a 3 + 3 a\ — 4 < 0, or 

• when 6 — 4 a 3 + a 4 < 0. 

In the former case Minkowski is the only vacuum solution and is stable, so that there is nothing further to worry 
about. However, in the latter situation it is possible to see, using the analogy (27), that the Minkowski vacuum is only 
a local minimum of an unbounded-from-below potential, and we should check its non-perturbative stability. Since c 
is constrained to be positive, a scenario of possible catastrophe is that a bubble of infinite negative energy forming 
inside of a Minkowski vacuum with c(r — > 0) — > +oo and c(r — > oo) — > 1. Repeating analysis similar to §IVD, we 
require §-s(c) > all the time, or with x = r^ 1 , 



where c(x) — c(r). Since c(x — > 0) — > 1, we can write 5=1 + (D(x a ) 7 a > 0, around x = 0, and then using s(l) = 1, 
(44) gives 



condition. 



F. Vacuum stability when there is only one classically stable vacuum 



d.r 



-s(c) > 



(44) 



dc 
xdx 




< 0. 



(45) 



That is, as x increases, c will keep decreasing, never being able to reach oo, so that the catastrophe will not occur. 
Therefore, a Minkowski vacuum, even when it seems to be a local minimum of an unbounded-from-below "potential" , 
is non-perturbatively stable. 
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V. HR BIGRAVITY 



We now consider the most general setup, where / is reinstated as a dynamical field, by having Mf finite. This also 
implies all five /3's are in play. 5 



A. Vacua 



As before, we work with the most general SO (4) -symmetric metric Ansatz in the Euclidean signature, 

g^dx^dx" = a g (r) 2 dr 2 + b g (r) 2 dQ 2 u , (46) 
f^dx^dx" = a f (r) 2 dr 2 + b f (r) 2 dn 2 u , (47) 

and use gauge freedom to fix a/ = 1. Each of (4) and (5) gives two independent equations, and there is one 
conservation equation (Bianchi identity). Barring redundancy, we end up with three equations for three yet-to-be 
determined functions, a g , b g and bf. 

^-( b ' g ) 2 + r ^^(^b 2 g + ^ 1 b g bf + 3^ 2 b 2 +^)a 2 g = °> ( 48 ) 

((3 1 b 2 g + 2p 2 b g b f + l3 3 b})(b g -a g b , f ) = 0. (50) 
From (50), again we have two separate branches of vacua. 



1. "CS" 

When ft b 2 g + 2 /3 2 b g b f + fa b 2 = 0, i.e., b g = X ±bf with ft X ± + 2 ft X± + h = 0, (48-49) are solved by 

, _ , , sin(V w (x±)r + 0i) 2 _ 2 l-«(x±)^ . . 

6 * =x±6 " bf = — 7m ' a ^ x± ^few (51) 

where 0i is an integration constant, and where we have defined 

m 2 M 2 

u(x) = -^iJ(/3 oa; 3 + 3A a ; 2 + 3/3 2 x + /3 3 ), (52) 
m 2 M 2 

v(x) = -^fW 1 x 3 + 3(3 2 x 2 + 3f3 3 x + f3 i ). (53) 

Note that both g- and /-metrics are (A)dS, with b>)R£ = 3 "_ ( f± } ^ and = 3v(x±)5£. We denote this branch 

x± 

as "CS" as it converges to the CS branch of dRGT in the limit Mf — » oo. 



2. "CM" 

When a g = b' g /b'f, the solutions for (48-49) are 



, - , , sm{<Jv(ci)r + 6 2 ) _ .... 

b g = Cl b f , bf = , a g = Cl , (54) 



§ Background solutions — in Lorentzian signature — we consider in this section were already found in [13, 14], and their perturbative analysis 
has been performed, although in different contexts, in [34, 37]. 
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where 6* 2 is an integration constant and c, is determined by 

v(ci) - u(ci) = 0. (55) 

In this case both metrics g^ v and describe spheres of radius Cj v(di)^ 1 ^ 2 and u(c"i) _1 / 2 respectively, that is (Anti)- 
de Sitter spaces of different curvature (analogous metrics were discussed in [37]). We denote this branch as "CM" 
because it converges to the CM branch in the limit Mf — > oo. We will discuss below how the two de Sitter metrics 
can be simultaneously written in cosmological coordinates with flat slicing and we will study perturbations on this 
background. 



B. de Sitter "CM" backgrounds and perturbations on top of them 

It is easy to see how we reach dRGT from HR by taking Mf — > oo limit. But obviously HR can have a much 
richer variety of vacua, including dS backgrounds in flat chart. To realize such backgrounds, we take (54) and choose 
9 2 = 7r/2. Also we constrain fj's such that v(di) > 0, and define 

v(c i ) = H 2 c 2 . (56) 

Then, (46-47), after a Wick rotation, becomes 

faAxW = cl ( - dt 2 + CO §ffi* d0 2 II ) , (57) 

f^dx" = -df'+ '^f d^, (58) 

which are conformal-to-dS metrics in the global coordinates. 

To analyze the perturbative consistency of this background, we use the planar coordinates for dS spacetime, and 
rescale t by i/cj to obtain 

fa = diag(-l,e 2Ht ,e 2Ht ,e 2Ht ), 



1 



U„ = ^dia^-l.e^e"",^'), (59) 

C i 

where Ci and H are determined by (55) and (56), i.e., 

-(A> + 3f +3| + |) = 0, (60) 



3^ 2 M 2 s( B 2 fa 

■'II ' m 2 s (q _ , na , , ; 



w2 ^( Ae , + 3 A + 3| + |) =0. (61) 

Note that, by going to conformal time r = — e~ Ht /H, it is easy to see that the metric g M „ describes de Sitter space 
with curvature H 2 , whereas / M „ gives de Sitter space with a different curvature (c"j if) 2 . 

We now expand the HR bigravity action, eq. (1), in metric perturbations. Referring the reader to Appendix A for 
algebraic details, the helicity-2 action is 

s hcl - 2 = / d 4 x {i(A TT ) 2 + \x^(K + 9 Jf)x^ + i(r TT ) 2 + ir TT (A - ™ 2 ff + ^)r TT } , (62) 

where X TT and Y TT are the linear combinations of the helicity-2 perturbations of g and /, A = e~ 2Ht A and where 

g2 J^2 _|_ 

m ' ff = cfWI 2 + M 2 ) c "' +2 ft c "' + W- ( 63 ) 
For the helicity-1 and -0 sectors, we obtain the following Hamiltonians 

Shel-! = / d 4 x {pV - i(p T ) 2 + i, T (A - m 2 ff + ^) q T } ■ (64) 

d 4 x|pg--p 2 + -( Z (A-m 2 ff + — Jqj, (65) 
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where all the non-physical gauge DOFs are solved away, and and are the surviving propagating helicity-O(-l) 
mode(s) and its conjugate momentum. 

To have a consistent theory, we need to make sure there is no ghost or tachyon. The absence of a tachyonic mode 
is guaranteed if m 2 s > 0. Inspecting the field redefinitions and canonical transformations performed for the metric 
perturbations, we see that under the ass umption of n o tachyon, the hazard of a ghost lurks only in the helicity-0 sector 
due to the appearance of the quantity \/m 2 s — 2H 2 . That is, the helicity-0 mode becomes a ghost when < 2 if 2 , 
which is the famous Higuchi ghost [38]. Therefore, requiring 

m 2 cS >2H 2 , (66) 

should be enough for HR bigravity on "CM" backgrounds to be consistent. 



C. An application: Vector modes in a more general "CS" branch 

Due to its gravity-modifying nature, it is natural to try to apply dRGT to cosmology. This, however, is not 
trivial, as flat cosmological solutions cannot be found for simple forms of the auxiliary metric [39]. Simple 
cosmological solutions can be found by looking for open spatial sections [11] (sec our discussion in section III A 
where these solutions appear naturally in the CS branch of our Euclidean construction). But such open solutions 
come with a rather worrisome feature that the vector perturbations on top of them do not propagate at the linear 
order [31], which is potentially dangerous because it could mean that they arc (infinitely) strongly coupled. Ref. [35] 
confirmed this concern by showing that once the background isometry is broken, an extra ghost mode appears. 
Additional confirmation was given in [32], which showed that the vector modes might get a kinetic term at higher 
order expansion. One may wonder if HR bigravity on a more general cosmological background can be free of such 
pathology. 

To check this possibility, we perform a double Wick rotation and a coordinate transformation on the metric of (51) 
with 6*i = 0, to obtain a cosmological background in HR bigravity: 

g„ v dx>*dx v = -dt 2 + ( sinh ^^-tY k 2 a ij dx i dx j , (67) 

v v M (x±) x± 1 

cosh 2 V" ( x±) ^ ,AI7?7Tx2 

f^dx^dx" = ^ — - dt 2 + ( sinh v } x± ' t) tfaadx'dx 3 , (68) 

Xl(l + 4§4sinh 2 V^MIA V«(X±) X± ' 

with <Tij = Sij — i+k™s 3 "x m x n ' x% ~ x > y> z - Since we are interested in the kinetic coefficients of the vector modes, let 
us focus on 

v = (4 Wi^) ' ^ = (4 v,cj + v,c t ) ' (69) 

where h and 8 are perturbations on top of the background defined by (67-68), Vj is a covariant derivative with respect 
to aij, and V 4 /iJj = = = V 4 £^ = V l ^ T . We can then perform a coordinate transformation to fix the gauge to 

ef = o. 

It is now straightforward to expand the equations of motion (4) for the metric <7 M „ in terms of (69), and at 0(h, 9) 
the Oi-component of (4) gives = 0. Using this result, the 0(h, #)-piece of the ^/'-component of (4) becomes 

v,c7 + v,C T = o, (70) 

i.e.,C7 = 0. 

In order to repeat the same algebra for the /-metric EOMs (5), it is convenient to perform another temporal 
coordinate transformation by 




||j sinh V^±)r + ^l + ^g|sinh 2 y^)r) , (71) 



13 



to get 

g^xW = C ° S ^ y^!_ -dr 2 + sinh y/^jr) V^W , (72) 

/ F d^d/ = -dr 2 + ( ^V^^ V^-cLrW , (73) 

v vMx±) ' 

so that perturbed equations of motion for the /-metric can be easily calculated. This transformation also affects the 
metric perturbations, and we will denote them as h^ v and 9^ v . Now, since (71) involves only the i-coordinate, 
and 0^ are different from h ^ and 0^ only by a factor of ^ , while $J and Q are the same as £j (that we have set to 
vanish by a spatial gauge transformation) and (J. That is, the transformation (71) does not mix up perturbations. 
Then after solving away 6^, the 0(h, #)-part of the zj-component of (5) turns out to be the same as (70). (To be 
precise, the full expression of the zj-component of (4) is (70) multiplied by a complicated function of t, whereas that 
of (5) is (70) multiplied by a slightly different function of r.) 

This result implies that, at the linear level, the equations of motion around the "CS" branch do not allow for 
propagating vectors, precisely as was found in the corresponding branch of the dRGT model. This is consistent with 
the results obtained in [34] . Since the dRGT vector modes obtain kinetic terms from the cubic expansion of the action, 
and dRGT is a special case of HR, it is straightforward to argue that the higher order expansion of the HR action 
will provide kinetic terms for its vector modes, so that even in HR the vector perturbations are strongly coupled on 
cosmological backgrounds. 



VI. CONCLUSIONS 



dRGT massive gravity has a rich and interesting vacuum structure. We have discussed the most general 50(4)- 
symmetric Euclidean vacuum solutions to the theory, recovering self-accelerating [11] as well Minkowski- like [22] 
metrics. 

After reviewing the spectrum of linearized perturbations about the Minkowski-likc backgrounds (those about curved 
backgrounds have been studied in detail, e.g., in [31, 32, 35, 40]), we have considered the possibility of non-perturbative 
transitions connecting different vacua. Remarkably, we have found that it is not possible to construct instantons 
supported by singular matter distributions if we require the equations of motion to make sense at least from a 
distributional point of view. Indeed, the continuity of the physical metric g^ v cannot be maintained without inducing 
curvature singularities in the fiducial metric while keeping a distributionally well-defined mass term. 

Then we moved on to discuss instanton transitions supported by regular matter distributions. Again, our result is 
negative: we have shown that there can be no smooth SO (A) instantons connecting the various vacua, unless we are 
willing to give up energy conservation or the null energy condition. More specifically, it is not possible to transition 
to/from a self-accelerating branch without violating energy conservation, while the null energy condition forbids any 
transition involving conformal-to-Minkowski backgrounds. 

The situation is different from the one discussed in [24] that analyzed instanton transitions in the DGP model. There 
it was shown that it is not possible to leave the self-accelerating branch without violating the null energy condition. 
Since the self-accelerating branch of the DGP model contains a ghost (that does violate energy conditions), it is then 
conceivable that a kink built out of such a ghost might provide the required negative energy [41]. In dRGT gravity 
a ghost is not sufficient to depart from the self accelerating branch, since, as we have seen, one needs an even more 
dramatic form of pathology: energy non-conservation. 

We have also discussed some aspects of the landscape of vacua in the Hassan-Rosen bigravity, constructing the full 
set of vacua that converge, in the limit where the auxiliary metric f^ v becomes non-dynamical, to the CS and CM 
solutions of dRGT discussed in section III A. In particular, we have shown that the branch converging to the CM 
background has well-behaved perturbations as long as the Higuchi bound is obeyed. We have also shown - consistently 
with the findings of [34] - that vector modes are non-dynamical in the branch of bigravity solutions that, in the limit 
Mf — > oo, converges to the self-accelerating CS branch. With this observation we can argue that this branch will 
feature strongly coupled (and possibly ghost-like) modes, in analogy with the situation discussed in [31]. 

For what concerns future work, it would be interesting to extend the instanton analysis to HR bigravity. Also, since 
it was shown in [35] that anisotropic backgrounds can generate ghost-like kinetic terms for some of the fluctuations 
about the CS solution, it might be possible to circumvent our no-go results of section IV by going beyond 50(4)- 
symmetric configurations. Even more important, of course, would be to find other vacua of these models with desirable 
features. 
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Appendix A: Perturbative expansion of HR model on the de Sitter "CM" background 

We expand the HR bigravity action, (1), in metric perturbations, 

where g and / are given by (59) and 

^oo = , h oi = hl L + e 2Ht d t B , hij = hf? + d^j + djg + e 2Ht {-25^ + 2d t d,E) , (A2) 
9 Q0 = -2cp, 6 0i = 6l + e 2Ht d t B , 6 l0 = eJ T + d l (J + d^J + e 2Ht (-25 tf w + 2didj£) , (A3) 

with dthl = dtfl = 0, Sijh?? = <5y0« T = °' Mi? = Ml? = and d £J = 9i@ = °- 

1. Helicity-2 

We define the linear combinations of hj? and 6j? by 

; tt — 2e m/2 ( r . x TT i M fy TT \ f) TT - 2eHt/2 ( xTT M q y tt\ , A4 n 

where the spatial indicies in A and Y are suppressed. Then A TT and Y TT do not mix in the helicity-2 part of (1), 
so that we obtain 

S hcl - 2 = / d 4 * {\{X^f + \X^(A + 9 Jf) X ^ + \(Y^f + (A - ml, + ^)y TT } , (A5) 

where m 2 s is defined in (63). Note that we have used (60-61) to simplify the result. 

2. Helicity-1 

We define new variables similar to the hclicity-2 case, 

,T C T _ ( - y-T , M f Y T \ „T VT _ ( X (U2) M g T \ 

J^Mjj + M) V M g ^'U Jc 2 M? + M 2 V c « M / V 

with the z-indicies suppressed in A^ 2 s, y,j 2 j, to split the helicity-1 part of (1) into X T and Y T sectors. 

In the X T sector, we introduce a gauge invariant variable vj^ by A^ = e Ht / 2 vJ^ + A^ — 2HXj 2 y to see that 
Xj 2 j completely disappears from the action, vj^ shows up as 

|d 4 ,{-^A(^ 1))2 }, (A7) 

i.e., = on shell, and there is nothing left in the A T sector. 
In the Y T sector, we first rescale Y T, s by Y^ 2 , = e lit / 2 uj 1 2 n to get 



'(2) 



5 ho i-i = J d A x Lhei-i = J d 4 x { - ^e 2H *U( r 1) (A - m^uj^ - jii^Aiif 

1 / Qff^ 1 \ 3 i 

-- (— - m? ff )wf 2) Awf 2) + «f 2 )^«a) ~ 2^ u a) Au (2)} • (A8) 
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Then we employ the Hamiltonian formulation and define the momentum P T conjugate to uj 2 ) by P T = 1 ■ 

Writing (A8) in J d i x(P T uJ 2 - ) — Hamiltonian) form, integrating out the non-dynamical constraint uj^ and performing 
a canonical transformation by 11 , 

rp y=Ar3P T /9P 2 2 \ T i T (-A)- 1 / 2 / T 3P T \ 

p =-n^{-Y p T + hr-"&> T }< "w = S^-( p + t 9 )' (A9) 

put the helicity-1 action in the canonical form: 

Shei-i = / d 4 x - \{p*Y + \<?(a m 2 cS + 9 | 2 ), T } . (A10) 

3. Helicity-0 

Again we combine the helicity-0 modes from the g-metric and those from the /-metric by 

<j>,B,1p,E = = =(ciX {1 2,3,4) + ^^(1,2,3,4)) , 

V% / ^"(1,2,3,4) M ff 



/ ^-(1,2,3,4) M s \ 

,/c?M 2 + M 2 V c * f 1 



(All) 



so that there is no mixing between X and y modes in the hclicity-0 action. 

All the X modes, just like the helicity-1 case, turn out to be either a gauge DOF or non-dynamical, so that none 
of them survives onshell. In the F-sector, we repeat the process of obtaining the canonical action for the helicity-1 
modes. That is, we first rescale Y's by 5^(1, 3,4) = e~ 3Ht / 2 cri,3,4 and F( 2 ) = e~ 7Ht / 2 a 2 to have 

Shd-o = J d 4 x I - 6H 2 a 2 - AHa l Ka 2 - ^TO 2 ff o- 2 Acr 2 + cri(4A - 6m 2 ff + 18H 2 )a 3 + 6Ha 2 Aa 3 
/ 27 \ 

-03 [2A - 6m 2 cS — J a 3 - 2(3H 2 - m^JffiA^ + (9P 2 - 4m 2 ff )er3Aer 4 (A12) 

-YlHo x a z - 4(T 3 Acr2 - 6o| + 4P(74Acri + 4cr 3 Acr 4 | . (A13) 

Then we introduce momenta, P3 and P4, conjugate to 0-3 and 04, and write (A12) in J d 4 x(P30- 3 + P4 04 — Hamiltonian) 
form: 

J d 4 x [p 3 <7 3 + Pi<?4 + vi{HP 3 + (4A - 6m 2 ff + 18H 2 )a 3 - 2(3P 2 - m 2 eS )Aa 4 } 
- lm 2 eS a 2 Aa 2 + a 2 (-e- 2Ht P 4 + 6HAa 3 ) - ^A" 1 P 4 - |p 4 A- 2 P 4 

-a 3 (2A - 6m 2 ff + ^-)<r 3 + (9P 2 - 4m^)<7 3 A <t 4 ] • (A14) 

(T2 is a non-dynamical constraint and can be integrated away, whereas o\ is a Lagrange multiplier whose elimination 
allows us to remove P 3 by P 3 = H- 1 { - (4A - 6™ 2 ff + 18P 2 )cr 3 + 2(3P 2 - m 2 ff )A<7 4 }. Thus we have only one 
dynamical conjugate pair (<t 4 , P 4 ), and with a canonical transformation by 



'hel-0 



2 A r m cSv /m 2 cS - 2H 2 



CT 4 



H I 2V3 

V3A- 1 (A m 2 ff t 3P 

m cS ^ml s - 2H 2 



q - (miff - 3H 2 )a 3 } , (A15) 



{*+(f-¥ + f fM- < A16 > 



^Here we assume rrig ff > 0, i.e., no tachyonic instability. If m^ s < 0, a slightly different canonical transformation is necessary, with a 
prefactor of «/— m'^ + A, which becomes imaginary at small enough scales. That is, at small enough scales the hclicity-1 modes have a 
negative definite Hamiltonian and become ghosts, just as seen in the example of dRGT perturbative analysis. 
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we finally get 



'hel-0 



9H< 



(A17) 
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